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A point particle of small mass m moves in free fall through a background vacuum spacetime metric 
g^ij and creates a first-order metric perturbation h^]^ that diverges at the particle. Elementary 
expressions are known for the singular m/r part of /iai.°* and for its tidal distortion determined by 
the Riemann tensor in a neighborhood of m. Subtracting this singular part h]^ from /i^^'' leaves 
a regular remainder /i^t ■ The self-force on the particle from its own gravitational field adjusts the 
world line at 0{m) to be a geodesic of g'^i, + h\^ . The generalization of this description to second- 
order perturbations is developed and results in a wave equation governing the second-order /i^5,° 
with a source that has an 0{mF') contribution from the stress-energy tensor of m added to a term 
quadratic in h]^^ . Second-order self- force analysis is similar to that at first order: The second-order 
singular field h?^ subtracted from h^^ yields the regular remainder h^ , and the second-order 
self- force is then revealed as geodesic motion of m in the metric g^^, + h + h^ . 

PACS numbers: 04.25.-g, 04.20.-q, 04.70.Bw, 04.30.Db 



I. OVERVIEW 

Recent, impressive fully relativistic numerical analysis 
has been brought to bear on a black hole binary system 
with a mass ratio of 100 to 1 [1, 2], and the evolution 
is followed for two full orbits before coalescence. The 
two disparate length scales of an extreme or interme- 
diate mass-ratio binary pose a challenge for numerical 
relativists to resolve the geometry in the vicinity of the 
small object while efhciently analyzing the remainder of 
spacetime and providing gravitational wave trains for a 
number of orbits. Second-order perturbation theory in 
general relativity might more efficiently meet the chal- 
lenge of the difficult numerical problems of extreme and 
intermediate mass-ratio binaries. 

Early descriptions of second-order perturbation the- 
ory [3-10] have focused on perturbations with no mat- 
ter sources and are typically limited to metrics with a 
substantial amount of symmetry. However, Habisohn 
[11] presents a fully general description of matter- free 
second-order perturbation theory for a background vac- 
uum spacetime metric g'^^^. 

Rosenthal [12-14] was first to describe a formal ap- 
proach to second order perturbation theory which in- 
cludes a small-mass (5-function point source. However, 
an actual application of his approach does not appear to 
be straightforward. 

The heart of this manuscript extends Habisohn's [11] 
second-order analysis to allow for a perturbing (5-function 
point mass. Our formalism is closely related to the tra- 
ditional description of linear perturbation theory. 

We begin in Section H with the formal expansion of the 
Einstein tensor, for a metric gab + ^a6, in powers of hab- 
First-order perturbation theory is summarized in Section 
HI for the case that the source is a (5-function object of 



small mass m. In the test mass limit m moves along a 
geodesic 70 of the background metric g^j,. With a finite 
mass m the metric is perturbed by the retarded field h}j^'^ 
at first order in m, and m's worldline deviates from 70 
by an amount of 0{m) as m itself interacts with h^l^^ as 
a consequence of the first-order gravitational self-force as 
described in Section IV. Throughout this manuscript we 
assume that the effects of m's spin and multipole struc- 
ture on its motion are insignificant when compared with 
the self-force effects. 

The extension of Habisohn's [11] second-order analy- 
sis to allow a (5-function point source demands careful 
consideration of the singular behavior of the metric in a 
neighborhood of m as described in Section V. Ultimately 
the wave equation for the second-order h"^]^^ appears in 
Eq. (26) as one might have expected, and the self-force 
analysis at second-order is seen to be similar in style to 
the analysis at first-order. 

The application of second order perturbation theory 
for a small mass still requires an effort which is strongly 
dependent upon the details of the actual problem of inter- 
est. Practical considerations are emphasized in Section 
VI. 



Notation and conventions 

In a neighborhood of a geodesic 70 of the background 
metric gJJ^ we use locally inertial and Cartesian (LIC) 
coordinates [15] where the timelike coordinate is t, the 
spatial indices i, j, k and / run from 1 to 3, the spatial 



coordinates are x* and r 



'77-y. 



In addition LIC 
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coordinates have special properties on 70: the coordi- 
nate t is the proper time, the spatial coordinates are all 
zero x^ — 0, the metric is the flat Minkowski metric r]ab, 
and all first coordinate derivatives of g'^f^ vanish. Second 
derivatives of g[]j on 70 determine a curvature length and 
time scale TZ, and the components of the Riemann tensor 
then scale as 1/7?.^ and their time derivatives along 70 



scale as 1/TZ^. After some fine-tuning of tlie coordinates 
[15-17], the metric in a neighborhood of 70 may be put 
into the form 

5°fc dx" dx^ = -qab da;" dx^ - x'x^R^t. (dt^ + 5ki dx^ dx^) 



where 



-xVi?°,^,dtdx'^+0(r■V7^3), 



(1) 



where the superscript " on the components of the Rie- 
mann tensor imphes that it is to be evaluated on 70. 
Also, both R%ij and R^^-jf are symmetric and tracefree 
in the indices i and j as consequences of the vacuum 
Einstein equations. 

Much of our analysis takes place in the buffer zone [16], 
a region spatially-surrounding 70 where m <^ r <^ TZ. In 
the buffer zone r is small enough compared to the cur- 
vature length scale, r ^ TZ, that the curvature of gJJ^ is 
barely apparent, and we have the luxury of being able to 
expand the actual metric g°^ -I- /ijj*^* away from flat space- 
time in powers of two simultaneously small numbers, m/r 
and r/TZ. 



II. EXPANSION OF THE EINSTEIN TENSOR 

We consider a perturbation hat of a given metric gab, 
and expand the Einstein tensor of the sum Gat [g + h) in 
terms of increasing powers of hah so that formally 

Gabig + h) = Gabig) + G^^ (g, h) + G^^^ {g,h) + --- (2) 

where Habisohn [11] describes an individual term in this 
expansion by 



Gi:^(.9,/^)-^ 



— -Gabig + \h) 



dA 



(3) 



A=0 



.in). 



This notation implies that the operator G^^^ (5, h) returns 
an expression that scales as (hab)"'- For n ~ 1 and gab 
being a vacuum solution of the Einstein equation, 

2Gi',^(.9, h) =. -V^V.hab - VaVhh\ + 2V^a'^^h), 

- 2Ra^b'-hcd + gabiWch^'d - W'^hcd), 

(4) 

where Vq is the derivative operator compatible with the 
metric gab- Habisohn [11] provides the following expres- 
sion for G\^{g, h) in his Eq. (3.1), 

+ (v[^/i'^l„) Vchdb - \C'' (2V(>fc)rf - Vdhab) 
+ [Ic^Cd - -^h^^V^V.Ka - \ {V^h-") ^ehcd 

+ ^h^''\/cCd + J (S/^h^^) \/chde] gab (5) 



Cd 



2V=/i, 



Vdhr 



(6) 



III. FIRST-ORDER PERTURBATION THEORY 
FOR A POINT MASS 

We next consider the consequences of adding an object 
of small size and small mass to, with to, ^ 7^, to the 
vacuum spacetime whose metric is 5°^. 

With a global coordinate system (T, X*), the stress- 
energy tensor for to- moving on a geodesic 70 of g^j, is 



Tabino) = TO —-6\X^ 



'PdT 



%iT)), 



(7) 



where 7o(T) gives the spatial position of the geodesic as 
a function of T, and the four- velocity Ua, a/— (7", and 
proper time r are all functions of T along the worldline. 
The dominant effect of Tab (70) on the spacetime metric 
results in the retarded metric perturbation h^j^^ propor- 
tional to m which solves 



Gabig" + h"""') = SnTabi-fo) + OK) 



(8) 



with appropriate boundary conditions. The superscript 
1 on any met 
for example. 



1 on any metric perturbation implies that /i^fc"* ^^ G(?n), 



^ab^ for an 0(to, ) metric 
perturbation and also use h'^ab = ^-'i^^ + ^"'^'^' 



Later we use 

. "" ^T + hlr + Oim^). 
For this linear perturbation problem, we expand the 
Einstein tensor in Eq. (8) using Eq. (2) and isolate the 
terms linear in m to obtain the first-order perturbation 
equation. 



:(i)^.o 



G'^>ig",h'n^S7rTabijo). 



(9) 



The Bianchi identity implies for arbitrary hab that if 
gab is a vacuum solution of the Einstein equation, then 



V"Gi^,'(.9,M 



0, 



(10) 



perhaps as a distribution. An integrability condition for 
Eq. (9) thus requires that ra6(7o) be divergence free. The 
assumption that the worldline of tti is a geodesic 70 of gj]^ 
guarantees that ^"'Tabi'fo) ~ and that the integrability 
condition is satisfied. 



IV. FIRST-ORDER GRAVITATIONAL 
SELF-FORCE 



After hl^f^*' is found using Eq. (9) there are several ways 
of calculating, understanding and interpreting the grav- 
itational self- force [17-23]. Our favorite is to note that 
^]iT^ is naturally decomposed within a neighborhood of 
7o into two complementary parts, 



Llrct _ ,is 
"-afc ~ ''■ab 



'^ab ■ 



(11) 



The first part h\f is the linear piece of the singular 
field ft,^j which is a special solution of 



Gabig"" + /i') = 8^r„6(7o) 



(12) 



with the notable features that /i^^: (1) may be expanded 
in powers of to, (2) is local to to and does not depend 
upon boundary conditions, (3) is accessible via an asymp- 
totic expansion [17-21] each term of which is singular or 
of limited differentiability on 70, and (4) does not exert a 
force on m itself, just as the Coulomb field of an electron 
at rest exerts no net force on the electron. 



The substitution ft,^^ = hlf + hlf + 0{m^) 



with hll 



0{m?), into Eq. (12) and the expansion of the Einstein 
tensor results in two equations, the first linear in m and 
the second quadratic. 






(13) 
(14) 



The inhomogeneous, linear singular field h^f looks like 
a Coulomb m/r field being tidally distorted by the Rie- 
mann tensor of g'^i^. We qualitatively describe h^f, using 
Lie coordinates associated with 70, as 



m 
r 



hlt^-{l + ^ 



X 

TV 



(15) 



only the scaling of the leading terms are shown, and this 
scaling is valid in the buffer zone, where to, <C r <C 7?.. We 
distinguish x from r to emphasize that x/r is generally 
finite but discontinuous C~^ in the limit r — > 0. The 
dominant term, scaling as just m/r, represents the linear 
in TO, terms in an m/r expansion of the Schwarzschild 
metric, as given in Eq. (A6) in Appendix A. The second 
term in the parentheses reflects the quadrupole distortion 
of the m/r field that is induced by the external Riemann 
tensor's tidal effects which scale as x^ /Ti? , as given by 
the terms proportional to to, in Eq. (A8). 

The complement of ft,^^ is the homogeneous regular 
field h]/^ = h]^^^ — h\f, from Eq. (11), which solves 



over a length scale comparable to to. An observer in a 
neighborhood of to,, with no a priori knowledge of the 
global spacetime, could measure the actual metric g^^ + 
h^ + h^^f at 0{m) and could distinguish the singular 
behavior of h^f from the remainder gj]^ + h^. However, 
the observer would be unable to distinguish h^ from 
(jfj^f, in the combination gj^^ + h]^ at linear order via local 
measurements only because gab~^^^ is a smooth solution 
of the vacuum Einstein equations at linear order. The 
observer would then naturally note that the worldline of 
TO is a geodesic 70 + 71R of the metric g°j + h^. The 
difference between the two worldlincs is denoted 71R and 
reflects the effects of what is often called the gravitational 
self-force, even though there is neither a force on m nor 
an acceleration of its worldline within the external metric 



.9°. 






It is apparent that an 0(771) coordinate transformation 
of the original LIC coordinates for 70 would remove the 
dipolc term in Eq. (17) and put the sum g'^f^ + h^ into 
the same form as displayed in Eq. (1), with 0{m) changes 
in the components of the external Riemann tensor. 

In an application h]/j^*' is typically found numerically 
while h^^ (or its approximation, cf. Section VI) is found 



analytically, then h]^ 



IR _ 



ab 



h\f gives the regular 



remainder (or its approximation) which is used to deter- 
mine the self-force and the appropriate geodesic 70 + 71R 

oig'ab 



LlR 
"afc ■ 



V. SECOND-ORDER PERTURBATION 
THEORY 

We assume that we have solved a first-order self-force 
problem of interest and have, in hand, ^;^h°*' ^af' ^^^ 
the initial geodesic 70 of gj]^ and the self-force modified 
geodesic 70 + 71R of g°j, + h'^J^. 

For the second order problem we also require /i^| which 
can be determined via an asymptotic expansion of Eq. 
(14), and scales as 



Gil^(3°,O=0. 



(16) 



The regular field h^ is smooth on 70 and, thus, qualita- 
tively described in a neighborhood of 70 by 



LlR 
'''ab 



m mx mx 



n n^ 7^3 



-f 



(17) 



with the LIC coordinates associated with 79. Each term 
takes the form of an external multipolc moment propor- 
tional to m. 

The regular field h]^ is added to gJJ^ to create the 
external metric 



cxt 
yab 



g'ab 



LlR 

'''ab 



(18) 



which governs the geodesic motion of to. After all, h^ 



l2S 
"-ah 



m 



X 

TV 



(19) 



with LIC coordinates. The dominant term, scaling as 
m? /r"^, is the term quadratic in to, in an m/r expan- 
sion of the Schwarzschild metric, as given in Eq. (A7) for 
n — 2. The second term in the parentheses reflects the 
quadrupole distortion of the m? /r"^ field that is induced 
by the external Riemann tensor's tidal effects which scale 
as x^ /TZ^ , as given by the 0{m'^) terms in Eq. (AS). 

To understand second-order perturbation theory re- 
quires understanding two distinct and critical roles 
played by the first order regular field h]/^. First, the 
stress-energy tensor of to is Tat (70 + 7ir), where the ar- 
gument implies that the worldline of to, is now a geodesic 



is a homogeneous solution of Eq. (16) with no variation of g^^^ -\- h^. The change in the stress-energy tensor re- 



suiting from the first-order self-force is 



Tab (70 +71r) - Tabi^o) 

= mA 



UaUbdr \^3[;^«_^.(y)] 



-gdT 



UgUb dr 



^iR^s'ix'-^oin (20) 



where the A operation reflects the 0(to) change in the 
quantity in parentheses which follows from changing the 



metric to gj^^ 



/^IR 
"ab 



from g^i^. Thus the difference be- 



tween the two stress-energy tensors is a distribution of 
0{m?) with support on 70 and consists of terms with a 
(5-function and with a gradient of a i5-function. 

A second effect of h^ on the second-order problem is 
the modification of the tidal environment of m by h^ 
which becomes an 0{m) part of the external metric as in 
Eq. (18). This creates 0{m) changes in the the external 
Riemann tensor's multipole moments. These changes are 
responsible for 0{m^) corrections to h^f which we label 
/i^j, . Thus the singular field is not derived solely from 
the initial geodesic and the background metric gJJ^, rather 
it specifically includes effects from the self-force modifi- 
cation of the geodesic and from the additional 0(?7i^) 
tidal distortion of h^^f caused by h^, and these 0{m^) 
contributions to the singular field constitute /i^jj, . 

The presence of h]^ in the external metric gj^^ + h]^ 
modifies the tidal effects of the external Riemann tensor 
on the singular field and Eq. (15) becomes 



An immediate application of this notation is in the 
recognition that 

GiV(5° + h'"", h'"" + h^""^) = SnTabijo + 7ir) + ^(m^), 

(23) 

which is the natural extension of Eq. (13) to second-order. 
The presence of h^ as part of the external metric in the 



.(1) 



first argument of G^^^ requires the addition of h^J to the 
second argument. We have already described h]^ in Eq. 
(11), and it is natural then to define h^ via 



i,2rct _ i,2R 
"ab — "-ab 



hlV 



l2S 
'^ab- 



(24) 



We now confront the second-order problem which re- 
quires a solution for h^]^*' from 

Gab{9° + h''"' + /i'"') = SnTabho + 7ir) + 0{m'), 

(25) 



when we are given the metric perturbations h^]^ , /t^jf, , 
h\f, h^^\ h'^f, and the worldhnes 70 and 70+71R. We ex- 
pand the left hand side about 5^]^, rearrange some terms, 

and substitute for G^^f, (5°, /i^'''^') from Eq. (9) to obtain 



:(l)/„0 L2rotN _ 



Gab (9^ h ) = ^^Tabho + 7ir) - 87rTab(7o) 



(26) 



as 



i2St 
^ab 



1 



7^2 



1 



m 

n 



where we are now using LIC coordinates for the geodesic 



7o + 7iR of g°j + h^^. The m/TZ term in the parentheses 



adds an 0{m^) contribution to /i^f,; however, the 0{m?) 
^2St jg naturally grouped with h\f because its presence 
in Eq. (21) algebraically resembles part of h\^ in Eq. (15) 
much more than any part of h^f in Eq. (19). 

Through second order the singular field is thus repre- 
sented by 



hlb 



'^ab 



hlV 



hit 



0{m^). 



(21) This wave equation for h^]^^ is the primary formal result 
of this manuscript. At the source each stress-energy term 
is 0{m); however, their difference is a distribution with 
support on 70 and is of 0{m^) as given in Eq. (20). 

The integrability condition for Eq. (26) is easily satis- 
fied away from 70 because there GJ^^ (g", /i^"^"') = and 
the fact that for any hab if G^j, (5°, h) = then it follows 

that V"G^j ((7°, h) = 0, as shown by Habisohn [11] in his 
Eq. (3.7). Thus the divergence of the right hand side is 
zero away from 70. The discussion of the integrability 
condition in a neighborhood of 70 is deferred until just 

(22) after Eq. (31) below. 



Eq. (26) becomes surprisingly transparent after some analysis (while cavalierly dropping terms of 0{rrv^) along the 
way) when h^^j^ is re-expressed with the substitutions /i;^^°* — h]^ + h]j^ and /i^^°' = h"^ + h^ + ^ab • Then the 
substitutions for the stress-energy tensors from Eqs. (13) and (23) lead to 

GiWg', h-^ + /.^st + /,2S) ^ GW(gO + /.!«, h'^ + h'^^) - Gi\\g", h^') d^{g\ h'^ + h'') . (27) 

Use of the identity in Eq. (B3) modifies the RHS with the result that 

I 



On the RHS, the fourth term cancels that part of the first 



term which is linear in h^^- 



The terms linear in h"^^ on 



the LHS and quadratic in h^^ on the RHS cancel from 
Eq. (14). The terms hnear in h"^^^ on each side of the 
equation cancel up to a term of 0{m^), which is ignored. 
When the dust has settled what remains is 



Gil^(.9°,0 = -Gil^(.9°>^^) 



(2)/'„0 l,lR^ 



(29) 



which reveals obvious consistency for this second or- 
der perturbation formahsm: When h]^, h^^ and K^ 
correctly capture their respective parts of the singular 
behavior of the retarded field, the regular remainder 
^^ab ~^ ^^ab app^ars as a source- free metric perturbation at 
second order in m as described by Habisohn [11]. The in- 
tegrability condition for Eq. (29) is satisfied in a manner 
similar to that for Eq. (26) away from 70. 

The second-order self-force is similar to the first-order 
self-force. In a neighborhood of m, h"^^^ is naturally de- 
composed into two complementary parts, /i^^°' = h^]^ + 
{h^^ + fi^f), where ft-^^ + h'^f exerts no force on m it- 
self. The second-order self-force then moves m along a 
geodesic of g|Jj, + h^f +^ hlf . 

The sanguine simplicity of Eq. (29) hides the complex- 
ity of its application. It might appear as though h"^ may 
be solved only in terms of h^ in a neighborhood of 70, 
but what is lacking is the description of the boundary 
condition which is typically given as a condition on the 
retarded field h^J^^ . To find h^ it is necessary first to find 
/i^™' and to evaluate h^f as an asymptotic expansion in 
a neighborhood of 70; these lead to h^^ — ft,;^^°* — h\f. 
With h^ the self- force modification of the worldline may 
be determined. At this point h^^^ and h^f are accessible 
via asymptotic expansions and /i^^°' could be evaluated 
via Eq. (26). Only then is h"^ able to be determined. 



VI. PRACTICAL CONCERNS 

In most situations, only an asymptotic approximation 
/i^f, to the exact ft,^^ is likely to be known, and as a conse- 
quence an actual application of the formalism described 
above is not as elementary as it might appear. In this 



case, Ky, 



/i^f, is an approximation to the actual 



regular field /i^j,. With these approximations some con- 
cerns appear in a neighborhood of the (5-function point 
source m. The proper evaluation of the self-force, via 
^ab' requires that /i^j, match both the value and first co- 
ordinate derivatives of h^i^ on 79. In turn, this requires 
that the difference ft-^^ — ft-^^ be zero on 70, and also, with 
Lie coordinates, that all first coordinate derivatives of 
this difference also be zero on 70. 

Experience [24-30] has shown that in numerical work 
if the difference ft-^^ — ft-^^ of these two singular fields is 
increasingly more differentiable, then the numerical anal- 
ysis will be increasingly more accurate. 

In some self-force analyses [31] 



"ab 



hll + 0{mx'^/rn'^) and 



We assume henceforth that we have such a precisely de- 
scribed approximation /i^j^ to /i^j,. 

For first order analyses, the integrability condition re- 
quired for using Eq. (9) to solve for h^]^^ is easily satis- 
fied. The approximation for ft,;^| is then accurate enough 
that h^^ is C^ on 70 , and the accuracy of the computed 
self-force effects are not limited by this approximation. 

To derive a second-order equation for /i^^ follow the 
same instructions as for Eq. (29) while using /ij^^ and ft-^^ 
instead of ft-J^^ and ft-^^, and do not use Eqs. (13), (14) or 
(23) for substitutions. The result is 

+ [SnTabho + 7ir) - GW(.g° + h'\ h'' + /i^^t)] 



[87rTabho)-G):>{g^\h'^)] 



W{^0 J,1SN 



(31) 



The integrability condition for using Eq. (31) to solve 
for h^^ is satisfied everywhere except, perhaps, precisely 
on 70 where the analysis entails some modest difficulty. 
The order terms associated with h}^l and hj^j + h^l (given 
above) provide an estimate for the behavior of the source 
on the righthand side in a neighborhood of 70. Most of 
the terms on the righthand side are either distributions 
or differentiable and well behaved on 70. The uncertainty 
involving the source is dominated by the 6)^^' (g°, /i^*^) and 

Gab (3°' ^'^'^^ + ^^'^) terms; each of these scales as two spa- 
tial derivatives of rri^x^ /r^TZ*^ which is 0{rn? x^ / r^Ti!^) 
and finite but discontinuous on 70. The divergence of 
this term is then 0{m?x/r^TZ'^) which diverges on 70. 
However, the integral of this divergence (contracted with 
a smooth test vector field of order unity) over a small 
volume of radius r* about m is then 0{m'^r'^/TZ'^). If we 
choose r* such that m, r*, and TZ are related by 



rl/n <m-^r^ 



<$:n 



(32) 



(30) 



then it follows that the integrated divergence over the 
volume of radius r^ is 0{m? /Ti?). For r > r^, the in- 
tegrability condition is satisfied. Thus, the integrability 
condition fails only at 0{m^) which does not hinder the 
analysis at O(to^). No fundamental difficulty prevents 
solving Eq. (31) for ft.^^. The resultant h^l is C^ on 70 
and is sufficient to find second order self-force effects. 



VII. SUMMARY AND CONCLUSIONS 

Upon reflection, Eq. (26) describes the second-order 
perturbation problem for a (5-function point mass in a 
quite satisfactory manner and is the primary result of 
this manuscript. The metric perturbation h'^]^^ may be 
determined directly, and the /i^j,, /i^j, decomposition of 
^ab* ^^ only required for determining the effects of the 
self- force. 

It is notable that the representation of a small mass m 
by a (5- function point source works as well at second-order 
as it does at first order. 
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Appendix A: Nonlinear perturbation theory and 
tidal distortion of a small black hole 

The simplest example of non-linear perturbation the- 
ory in General Relativity involves perturbing flat space- 
time by putting a small, spherical object of mass m down 
on the origin of Minkowski space. Outside the object 
the geometry must be the Schwarzschild metric from 
Birkhoff 's theorem. 

The usual coordinates of Minkowski space form an LIC 
coordinate system because the spatial origin x* = is a 
geodesic, and the other LIC conditions are clearly satis- 
fied. We define a covariant vector in the radial direction 
via Ui — ViT. With a Schwarzschild black hole of mass 
m present at the spatial origin, the metric takes the un- 
familiar form 



9ab 



^''"dx^dx'' = - f 1 - ::::: 



2m 
r 



dt' 



2m 



nkUidx dx 



+ {Ski - nkni)dx^ dx'- . 



An alternative description of this 
Schwarzschild metric is 



form 



schw 

9ab - Vab 






(Al) 
of the 

(A2) 



where o^a6 is to be identified as the singular field from 
self-force analysis, and the leading subscript implies 
that this monopole part of the singular field is spherically 
symmetric. From Eq. (Al) it follows that 



o/i^bdx"dx* 



/ schw 
\9ab 



2m 



-dt' 



riab)dx''dx'' 
2m 
r — 2m 



Ukfiidx dx 



The nth order part of o^a6 scales as m^' 
isolated with 



oK 



nS 

6 = 



— — /)S 
dm" "^ 



(A3) 
and may be 

(A4) 



m=0 



This provides the formal representation 

oo 

hlb = E ''ab ■ (A5) 

n=l 

For our elementary example, the first term in this sum is 

ohlfdx'^dx'' = —dt^ + —nknidx^dx\ (A6) 

r r 



oh^h dx"dx 



— ] nknidx''dx'-. (A7) 

r ' 



In this treatment of the Schwarzschild metric the singu- 
lar features of oh^f are identified, and the absence of a 
regular field /i^j, is assured by the flat nature of the initial 
Minkowski metric. 

A more subtle example places a Schwarzschild black 
hole in a region of spacetime that is empty but has slowly 
changing curvature from some distant source. In that 
case the metric of a black hole placed on the origin of the 
LIC coordinate system of Eq. (1) would be perturbed 
by the background curvature and could be analyzed by 
use of the Regge- Wheeler [32] formalism. The boundary 
condition at large r requires that the perturbed metric 
approach the form given in Eq. (1). The boundary condi- 
tion as r — > 2m requires that the perturbation be well be- 
haved on the future event horizon of the small black hole. 
In the time independent limit the wave equations for the 
metric perturbations admit analytic solutions which sat- 
isfy the boundary conditions [17]. 

The dominant tidal effects present in both h]^ and h?^ 
are seen in the quadrupole I — 2 terms of Eq. (9) of [17], 
which we reproduce here as 



2"q6 



hl^.dx'^dx^ ^ Rl^tjx'x^[{4:m/r-Am^/r^)dt^ 
+ 2m^/r^{5ki - nkUi) dx dx ] 

0{mx^/rn^ 



+0{m^x^/rn^) + 0{m^x^/r^n^) (A8) 



The order terms here result from the possible slow time 
dependence of the tidal field and are all much smaller in 
the buffer zone than the explicit terms provided. 

A more extensive analysis of ft,^^ in a similar style is 
given in [20] . An alternative treatment in a dramatically 
different style is given in [21]. 



Appendix B: Useful Identity 

An identity used in deriving Eqs. (29) and (31) results 
from considering two different expansions of the same 
expression G{g^ + h^^ + h^^). On the one hand, treating 
h^ + h^f as a single quantity, it expands to be 

GiWg', h'"" + h'') + G^^h9\ h'"" + h'') + 0{m^). 

(Bl) 

On the other hand, first grouping h]i^ with g^j^ while 
expanding in powers of h]^ , and subsequently expanding 
in powers of h]^ , it becomes 

G«(5°, /.i^) + G^^{g\ h'^) + d^{g\ h'^) 

+ Gi6'(5° + /^'^/^"^) + 0(m3). (B2) 



Equating these two expressions reveals that 

+ Oim^). (B3) 
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